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Abstract 

The high-temperature expansion of the spin-spin correlation function of the two-dimensional 
classical XY (planar rotator) model on the square lattice is extended by three terms, from order 
21 through order 24, and analyzed to improve the estimates of the critical parameters. 
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Tests of increasing accuracy- of the BKT theory^ of the two-dimensional XY model 
critical behavior have been made possible by the steady improvements of the computers 
performances and the progress in the numerical approximation algorithms. However, the 
critical parameters of this model have not yet been determined with a precision comparable 
to that reached for the usual power-law critical phenomena, due to the complicated and 
peculiar nature of the critical singularities. Therefore any effort at improving the accuracy of 
the available numerical methods by stretching them towards their (present) limits should be 
welcome. After extending the high-temperature(HT) expansions of the model in successive 
steps^ from order (3^^ to jS"^^, we present here a further extension by three orders for the 
expansions of the spin-spin correlation on the square lattice and perform a first brief analysis 
of our data for the susceptibility and the second-moment correlation-length. More results and 
further extensions both for the square and the triangular lattice^ will be presented elsewhere. 
Our study strengthens the support of the main results of the BKT theory already coming 
from the analysis of shorter series and suggests a closer agreement with recent high-precision 
simulation studies^i^ of the model. 

The Hamiltonian 

H{v} = -2 J J2 ^(^ ■ ^(^) (1) 

nn 

with v{r) a two-component unit vector at the site r of a square lattice, describes a system 
of XY spins with nearest-neighbor interactions. 
Computing the spin-spin correlation function, 

C(0,x;/3) =< s(0) ■ >, (2) 

(for all values of x for which the HT expansion coefficients are non-trivial within the max- 
imum order reached), as series expansion in the variable (3 = J/kT, enables us to evaluate 
the expansions of the l-th order spherical moments of the correlation function: 

m(')(/3) = ^|f|'<s(0)-s(x)> (3) 

X 

and in particular the reduced ferromagnetic susceptibility xiP) = m^^^P). In terms of 
m^'^\(3) and x{(^) we can form the second-moment correlation length: 

e^(/3)=m(2)(/?)/4x(/?). (4) 
Our results for the nearest-neighbor correlation function (or energy E per link) are: 



^ 2^ 3 48 120 1440 5040 645120 

6506950039 17 1102473407093 6986191770643 21 

26127360 2612736000 ^ 14370048000 

1657033646428733 

4138573824000 ^ ^ 



(5) 
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For the susceptibility we have: 

X = 1 + 4/5 + 12/32 + 34/33 + 88(3' + + 529/36 + ^^(3' + + ^^^p^ 

3 12 2 oU 

2608499^10 3834323 ^„ 1254799^12 84375807^1, 6511729891^14 

H B^^ H B^^ H B^^ H B^-^ H B^' 

180 120 ^ 18 ^ 560 ^ 20160 ^ 

66498259799^1, 1054178743699 ^i^ 39863505993331 17 

+ -\ B -\ B 

96768 ^ 725760 ^ 13063680 ^ 

19830277603399 18 8656980509809027 2985467351081077 20 
^ 3110400 ^ 653184000 ^ 108864000 

811927408684296587 21 399888050180302157 22 245277792666205990697 23 
^ 14370048000 ^ ^ 3448811520 ^ ^ 1034643456000 ^ 



83292382577873288741 ^ , 

+ /3 +0f/3 ) 

172440576000 m ^ ; 



(6) 



For the second moment of the correlation function we have: 

m2 = 4/3 + 32/32 ^ ig2/33 + 672/3^ + ^/3^ + '^p' + ^^/3^ + 77996/3« 

13484753^0 28201211 ^ 611969977 202640986^12 58900571047^1, 

H /3^ H 0^ H /3^^ H /3^2 _^ gi3 

60 45 360 45 5040 

3336209179^14 1721567587879^1, 16763079262169 ^i^ 5893118865913171^17 

H /3^^ H /3^^ H /3^6 H 

112 ^ 23040 90720 13063680 

17775777329026559 ^10 1697692411053976387 ^iq 41816028466101527 

H B^^ H /3^^ H /320 

16329600 653184000 ^ 6804000 ^ 

206973837048951639371 21 721617681295019782781 22 
^ 14370048000 ^ 21555072000 ^ 

79897272060888843617033^23 2287397511857949924319^24 ^/^25n 
H /323 H B^' + 0(3^^) (7) 

1034643456000 ^ 12933043200 ^ ^'^ ^ ^ ^ 

The coefficients of order less than 22 were already tabulated in Refs.-, but for complete- 
ness we report all known terms. As implied by eq.([T]), the normalization of these series 
reduces to that of our earlier papers^ by the change /3 —>■ (3/2. 

Let us now list briefly the main predictions^, of the BKT renormalization-group analysis 
to which the HT series should be confronted in order to extract the critical parameters. 

As /3 — /3c, the correlation length ^'^{(3) = m^'^\(3) / Ax{(3) is expected to diverge with the 
characteristic singularity 

e'(/5)(xd(/3) = ea;p(6/O[l + 0(r)] (8) 

where t = 1 — (3/ (3c- The exponent a takes the universal value a = 1/2, whereas 6 is a 
nonuniversal positive constant. At the critical inverse temperature (3 = (3c, the asymptotic 
behavior of the two-spin correlation function as = r ^ 00 is expected^ to be 

<.(0)..(f)>oc(!^|H- 0(5^^)1 (9) 
Universal values i] = 1/4 and = 1/16 are predicted also for the correlation exponents. 
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A simple non-rigorous argument based on eqs. ([8]) and ([9]) suggests that, for / > i] — 2, 
the spherical correlation moment m^^\(3) diverges as r ^ 0^ with the singularity 

mW(/?) (X r-^C2-,+^(/3)[l + Oir^'Hnr)] (10) 

This argument was challenged^ by a recent renormalization group analysis implying that 
the logarithmic factor in eq.([9]) gives rise to a less singular correction in the correlation 
moments, taking, for example in the case of the susceptibility, the form 

m^'\(3)^ea;\m + cQ] (11) 

where Q = 2(in{l)+u)'2 + ^i^^iO'^^) m is a non universal parameter. 
By eqs.® and (ITOl) . the ratios rn{m''^^) = anV'^i+i of the successive HT expansion 
coefficients of the correlation moment m*^'^(/5), for large n should behave^ as 

r„(m«) =P, + Q/{n + 1)^ + 0{l/n) (12) 

with C = 1/(1 + cr), to be contrasted with the value C = 1 which is found for the usual 
power-law critical singularities. 

To begin with, let us assume that a = 1/2 as expected, so that ( = 2/3. FigH] gives 
a suggestive visual test of the asymptotic behavior of some ratio sequences rn{m^''^) by 
comparing them with eq. (fT2!) . The four lowest continuous curves interpolating the data 
points are obtained by separate three-parameter fits of the ratio sequences r„(x), r„(m'^^/^)), 
rn{m^^^) and r„(m(^)) to the asymptotic form a + b/(n + 1^^^ + c/{n + 1) of eq. (fT2l) . In the 
same figure, the two upper sets of points are obtained by extrapolating the alternate-ratio 
sequence for the susceptibility, first in terms of l/(n + 1)^/^ and then in terms of l/(n -f- 1). 
The values of jSc indicated by the fits of the ratio sequences, range between 0.5592 and 
0.5611. 

A more accurate analysis can be based on the simple remark that, near the critical point, 
by eq.([8]) and ffTOj) (or eq. ffTTj) ). one has ln(x) = ci/r'^ + C2 + ... Therefore, if a = 1/2, 
the relative strength of the I/a/t and 1/r singularities in the function L{a, (3) = (a + 
ln(x))^ is determined by the value of the constant a. If we choose a ^ 1.19, the function 
L{a, (3) is approximately dominated by a simple pole and we can expect that the differential 
approximants (DAs)^ will be able to determine with higher accuracy not only the position, 
but also the exponent of the critical singularity. Using inhomogeneous second-order DAs 
of L{a,j3), we can locate the critical singularity at f3c = 0.5598(10). By analysing in the 
same way the series data truncated to order 21 which were previously available , we would 
get the estimate Pc = 0.5588(15). A consistent estimate j3c = 0.558(2) had been obtained 
in earlier independent^- studies of the same series using Fade approximants or first-order 
DAs. Older studies^ of slightly shorter series also indicated values of j3c in the same range, 
but with notably larger uncertainty. Thus our new series results indicate a stabilization 
and a sizable reduction of the spread for the (3c estimates. Our uncertainty estimates are 
generally taken as the width of the distribution of the values of (3c in the appropriate class of 
DAs. FigE] shows the singularity distribution (open histogram) of the set of quasi-diagonal 
DAs which yield our new estimate. These are chosen as the approximants [k,l,m;n\ with 
17 < k + l + m + n < 22. Moreover, we have taken — /|, \l — m\, \ k — m\ < 3 with fc, /, m > 3 
and 1 < n < 7. The class of DAs can be varied with no significant variation of the final 
estimates, for example by further restricting the extent of off-diagonality, or by varying the 
minimal degree of the polynomial coefficients in the DAs. No limitations have been imposed 
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on the exponents of the singular terms or on the background terms in the DAs in order to 
avoid biasing the Pc estimates. Should we require that the exponent of the most singular 
term in the approximants differs from -1, for example, by less than 20%, we would obtain 
/3c = 0.5602(5), well within the uncertainty of our previous unrestricted estimate. The 
vertical dashed line in FigI2] shows the value Pc = 0.55995 suggested by the simulation of 
Ref.-. Although no explicit indication of an uncertainty comes with this estimate, an upper 
bound to its error might be guessed from the statement^ that the simulation can exclude 
values larger or equal than jSc = 0.56045 for the inverse critical temperature. 

Biasing with Pc = 0.5598(10) the set previously specified of second-order DAs of L{a,P), 
leads to the exponent estimate a = 0.50(1). Figj2]also shows the distribution of the exponent 
estimates (hatched histogram) from this biased set. The uncertainty we have reported for a 
accounts not only for the width of its distribution shown in Figl2|, but also for the variation 
of its central value as the bias value of Pc is varied in the uncertainty interval of the critical 
inverse temperature. Essentially the same value of a would be obtained from the analysis 
of a series truncated to order 21. 

While, as one should expect, the DA estimate of Pc is rather insensitive to the choice of 
a, the estimate of the exponent a and the width of its distribution are fairly improved by 
our choice of a. Taking for example a = 0, we would find a = 0.53(4), which shows how the 
convergence of the exponent estimates is slowed down by the more complicated singularity 
structure of L{0,P). Similar values of a were found in previous studies of shorter series. 
Probably for the same reason, also the central values of the t] estimates obtained from the 
usual indicators are still slightly larger than expected. For example, by studying the function 
H{P) = ln(l + m^^''/x^)/ln(x) (or analogous functions of different moments), we can infer 
rj = 0.260(10). The function'' D{P) = \n{x) — (2 — ri)\'n.{^) and its first derivative are also 
interesting indicators of the value of rj. Taking rj = 1/4, Fade approximants and DAs do 
not detect any singular behavior of D{P) or of its derivative a.s P —>■ Pc, thus confirming the 
complete cancellation of the leading singularity in D{P). Moreover, this behavior seems to 
exclude the form eq. fllOl) of the corrections which implies the presence of weak subleading 
singularities, while it is compatible with eq. ffTT]) . 

In conclusion, our analysis suggests that, in spite of their diversity, the HT extended series 
approach and the latest most extensive simulation are competitive and lead to consistent 
numerical estimates of the highest accuracy so far possible. 
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FIG. 1: Ratios of the successive HT-expansion coefficients vs. l/(n + 1)^/^: for the susceptibility 
X (open circles), for m*^^/^-* (rhombs), for m^^^ (squares) and for m*-^^ (triangles). The four low- 
est continuous curves are obtained by separate three-parameter fits of each ratio sequence to its 
leading asymptotic behavior eq. (jl2p . The data points represented by crossed circles are obtained 
by extrapolating the sequence of the susceptibility alternate ratios with respect to 1/ n^/^ , and the 
continuous line interpolating them is the result of a two-parameter fit of the last few points to 
the expected asymptotic form a + b/n. The small black circles are obtained by a further extrap- 
olation of the latter quantities with respect to 1/n. The continuous line interpolating the black 
circles is drawn only as a guide to the eye. The horizontal broken line indicates the critical value 
Pc = 0.55995 suggested by the simulation of Ref.— 
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FIG. 2: Distribution of singularities for a class of second-order inhomogeneous DAs of L(1.19, /3) = 
(1.19 + Inx)^ versus their position on the P axis(open histogram). The central value of the open 
histogram is Pc = 0.5598(10). The bin width is 0.0007. The vertical dashed line shows the critical 
value Pc = 0.55995 indicated by the simulation of Ref.'' for which one can guess an uncertainty at 
least twice smaller than ours. The hatched histogram represents the distribution of the exponent 
a obtained from DAs of L(1.19, /3) biased with Pc = 0.5598, vs. their position on the a axis. The 
central value of the hatched histogram is o" = 0.500(1) and the bin width is 0.0015. The variation 
of the central value of a as Pc varies in its uncertainty interval is 0.01. This value can be taken as 
a more reliable estimate of the uncertainty of a. 
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